Phase fluctuations in atomic Bose gases 
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We improve on the Popov theory for partially Bose-Einstein condensed atomic gases by treating 
the phase fluctuations exactly. As a result, the theory becomes valid in arbitrary dimensions and is 
able to describe the low-temperature crossover between three, two and one-dimensional Bose gases, 
which is currently being explored experimentally. We consider both homogeneous and trapped Bose 
gases. 
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Introduction. — One of the most important features 
of the trapped Bose-Einstein condensed atomic gases is 
the remarkable control that can be achieved experimen- 
tally over the relevant physical parameters. As a result 
these gases can be studied in various physically distinct 
regimes, where their properties are quite different. The 
latest achievement in this respect is the experiment by 
Gorlitz et al. Q , in which one and two-dimensional Bose- 
Einstein condensates were created by reducing the tem- 
perature and the average interaction energy of the atoms 
below the energy splitting of either two or one of the di- 
rections of the harmonic trapping potential, respectively. 

As is well known, the physics of one and two- 
dimensional systems is fundamentally different from 
the physics in three dimensions. This difference is 
for instance illustrated by the famous Mermin-Wagner- 
Hohenberg theorem |^,|| , which states that in one dimen- 
sion Bose-Einstein condensation in a homogeneous Bose 
gas never occurs, whereas in two dimensions a conden- 
sate can only exist at zero temperature. In both cases, 
this is due to the enhanced importance of phase fluctua- 
tions. The Mermin-Wagner-Hohenberg theorem is valid 
only in the thermodynamic limit and does not apply to 
finite-size systems. In one or two-dimensional trapped 
Bose gases a condensate can therefore exist if the exter- 
nal trapping potential sufficiently restricts the size of the 
atomic gas cloud |^,^. 

The physics of low-dimensional Bose gases is in fact 
even more interesting, because, notwithstanding the 
Mermin-Wagner-Hohenberg theorem, a dilute homoge- 
neous two-dimensional Bose gas is expected to undergo 
at a nonzero critical temperature a true thermody- 
namic phase transition, that is known as the Kosterlitz- 
Thouless transition Q. Below the critical temperature, 
the gas is superfluid but has only algebraic long-range or- 
der. This so-called topological phase transition is there- 
fore not characterized by a local order parameter, but 
by the unbinding of vortex pairs and the resulting de- 
struction of superfluidity. Since the Mermin-Wagner- 
Hohenberg theorem forbids a true Bose-Einstein conden- 
sate in two dimensions, the superfluid phase is only char- 
acterized by the existence of a so-called "quasiconden- 



sate" . This important concept was first introduced by 
Popov and roughly speaking corresponds to a con- 
densate with a fluctuating phase. 

Although introduced theoretically as early as in 1972, 
the actual observation of such a quasicondensate has 
only very recently been made in a spin-polarized atomic 
hydrogen adsorbed on a superfluid ^He surface by Sa- 
fonov et al. In particular, this experiment measures 
the three-body (dipolar) recombination rate of a spin- 
polarized atomic hydrogen gas and observes a reduction 
in the associated rate constant due to the presence of a 
quasicondensate. Qualitatively, this reduction was an- 
ticipated by Kagan, Svistunov and Shlyapnikov, but the 
magnitude of the effect turns out to be much larger than 
predicted jo) . The reason for this discrepancy is still not 
fully understood, although a physical mechanism for the 
additional reduction was already suggested by Stoof and 
Bijlsma before its observation [p^ . 

In the context of trapped atomic Bose gases, the pos- 
sibility of observing a quasicondensate has been explored 
by Petrov et al. |ll| , p^ . Although this presents an im- 
portant first step towards understanding the physics of 
low-dimensional Bose gases, the approach can only be 
justified close to zero temperature. The reason for this 
is that density fluctuations are neglected from the out- 
set. As a result, depletion of the quasicondensate, due to 
either quantum or thermal fluctuations, cannot be prop- 
erly accounted for. In particular, the approach does not 
lead to an equation of state for the Bose gas. 

The main aim of the present Letter is to overcome this 
problem and to formulate a microscopic theory that in- 
cludes both density and phase fluctuations of the Bose 
gas. It is similar in spirit to the successful Popov theory 
for three-dimensional Bose-Einstein condensed gases, and 
can be used to study in detail the thermodynamic behav- 
ior of low-dimensional degenerate Bose gases. Moreover, 
it describes also the crossover between three, two and 
one-dimensional Bose gases. To explain most clearly how 
this can be achieved, we first discuss the homogeneous 
case. After that, we generalize the theory to inhomoge- 
neous gases. 

Modified Popov theory. — Formulating a microscopic 
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theory for lower-dimensional Bose gases is complicated 
by the fact that mean-field theory is plagued with in- 
frared divergences. To facilitate the discussion of how 
to deal with these divergences, we first recapitulate the 
expressions for the density n and the chemical potential 
/X that follow from the usual Popov theory for partially 
Bose-Einstein condensed gases. For a Bose gas in a box 
with volume V they read |^ 13 
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Here uq is the density of the Bose-Einstein conden- 
sate, Vb^(x — x') is the bare two-body interaction po- 
tential, ek = ?i^k^/2m is the kinetic energy of the atoms, 
Tiwk = (ck + 2noVbek)^^^ is the Bogoliubov dispersion 
relation, and N{x) = l/{e/^ — 1) is the Bose-Einstein 
distribution function where (3 = l/k^T is the inverse 
thermal energy. 

In agreement with the Mermin-Wagner-Hohenberg 
theorem, the momentum sums in Eqs. (0) and con- 
tain terms that are infrared divergent at all tempera- 
tures in one dimension and at any nonzero temperature 
in two dimensions. The physical reason for these "dan- 
gerous" terms is that the above expressions have been de- 
rived by taking into account only quadratic fluctuations 
around the classical result ng, i.e., by writing the annihi- 
lation operator for the atoms as '0(x) = -y«o + 
and neglecting in the hamiltonian terms of third and 
fourth order in ip'{x.). As a result the phase fluctua- 
tions of the condensate give the quadratic contribution 
no(x(x)x(x)) to the right-hand side of the above equa- 
tions, whereas an exact approach that sums up all the 
higher-order terms in the expansion would clearly give 
no contribution at all to these local quantities because 
(^^-^x^^) ^m^)^ = 1 + (x(x)x(x)) + . . . = 1. To correct for 
this we thus need to subtract the quadratic contribution 
of the phase fluctuations, which from Eqs. (|l|) and (|^) is 
seen to be given by 
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As expected, the infrared divergences that occur in the 
one and two-dimensional case are removed by performing 
this subtraction. 

After having removed the spurious contributions from 
the phase fluctuations of the condensate, the resulting 



expressions turn out to be ultraviolet divergent. These 
divergences can be removed by the standard renormal- 
ization of the bare coupling constant Vq. Apart from a 
subtraction, this essentially amounts to replacing every- 
where the bare two-body potential Vq by the two-body 
T-matrix evaluated at zero initial and final relative mo- 
menta and at the energy — 2/z, which we denote from now 
on by T^^{—2fi). Note that the energy argument of the 
T-matrix is ~2fi, because this is precisely the energy it 
costs to excite two atoms from the condensate 
this manner, we finally arrive at 
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/i = (2n - no)r^^ (-2/i) = {2n' + m)T''' (-2/i) , (5) 

where rt' = ti — rig represents the depletion of the conden- 
sate due to quantum and thermal fluctuations and the 
Bogoliubov quasiparticle dispersion now equals ?iWk = 

1 /2 

[e^ -I- 2noT^^(— 2/x)ek] . The most important feature 
of Eqs. (|^) and (^) is that they contain no infrared and ul- 
traviolet divergences and therefore can be applied in any 
dimension and at all temperatures, even if no condensate 
exists. How this can be reconciled with the Mermin- 
Wagner-Hohenberg theorem is discussed next. 

One dimension. — To understand the physical mean- 
ing of the quantity no in Eqs. (§ and (|), we must de- 
termine the off-diagonal long-range behavior of the one- 
particle density matrix. Because this is a nonlocal prop- 
erty of the Bose gas, the phase fluctuations contribute 
and in the limit |x| cx), we find 
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Moreover, using Eq. and carrying out the renormal- 
ization of the bare coupling constant, we obtain 
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At zero temperature, the quantity ([x(x) — x(0)] ) di- 
verges logarithmically for large distances, which leads to 
algebraic off-diagonal long-range order in the one-particle 
density matrix. In detail we can show that the leading be- 
havior of the zero-temperature one-particle density ma- 
trix for |x| ^ ^ is 



(V't(x)v;(o)) 



no 



(8) 



where rj = l/47rno^ is the correlation- function exponent 
and ^ = ?i/[4ninoT^^(— 2^)]^/^ is the correlation length. 
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At nonzero temperatures, ([x(x) — x(0)] ) diverges lin- 
early for large distances and the one-particle density ma- 
trix thus no longer displays off-diagonal long-range order. 

A few remarks are in order at this point. Most impor- 
tantly for our purposes, the asymptotic behavior of the 
one-particle density matrix at zero temperature proves 
that the gas is not Bose-Einstein condensed and that 
no should be identified with the quasicondensate den- 
sity. Moreover, from our equation of state we can show 
that in the weakly-interacting limit Aim^ ^ 1, the frac- 
tional depletion of the quasicondensate is {n — no)/n — 
(7r\/2/4— l)/47rn^ and therefore very small. Keeping this 
in mind, Eq. is in complete agreement with the exact 
result obtained by Haldane |T^. Note that our theory 
cannot describe the strongly- interacting case Attti^ ^ 1, 
where the one-dimensional Bose gas behaves as a Tonks 
gas ll7[|l|]. Finally, our results show that at nonzero 
temperatures not even a quasicondensate exists and we 
have to use the equation of state for the normal state 
n = N{ei^ + hT, — ^)/V to describe the gas. Here, the 
Hartree-Fock self-energy satisfies HT, = 2nT^^(— 

Two dimensions. — Applying the same arguments 
in two dimensions leads to the conclusion that at 
zero temperature, uq corresponds to the condensate 
density, whereas at a nonzero temperature, it repre- 
sents the quasicondensate density. In particular, the 
correlation-function exponent is 77 = 1/noA^ where A = 
(27r?i^/mfcBT)^/^ is the thermal de Broglie wave length. 
Due to the mean-field nature of the modified Popov the- 
ory, the Kosterlitz-Thouless transition is absent and a 
nontrivial solution of the equation of state exists even if 
r] > 1/4. This can be corrected for by explicitly including 
the effect of vortex pairs in the phase fluctuations. As 
we show in a future paper, this is achieved by using the 
modified Popov theory to determine the initial values of a 
renormalization-group calculation for the superfluid den- 
sity and the fugacity of the vortices. It should, however, 
be noted that for many applications we are not interested 
in the phase fluctuations and this additional renormaliza- 
tion is not very important. This is for instance true for 
the reduction of the three-body recombination rate con- 
stant of a hydrogen gas ||l9|] . 

At zero temperature, the fractional depletion of the 
condensate in the Popov approximation was first calcu- 
lated by Schick @ and is T'^^{~2^) /An where the chem- 
ical potential satisfies /i — nT^^(— 2/i). The correspond- 
ing result based on Eq. (0) is (1 - In 2) r2^(-2^)/47r 
where ^ now satisfies Eq. (0). Thus the depletion is re- 
duced by a factor of approximately 3. In two dimensions, 
the two-body T-matrix cannot be approximated by an 
energy-independent constant and this is the reason why 
we have been very careful about the appropriate energy 
of the two-body collisions pQ] . The T-matrix now obeys 



where a is the scattering length. 

Three dimensions. — We know that the Popov the- 
ory has been very successful in describing the properties 
of three-dimensional trapped Bose gases. It is therefore 
important to mention that, although the modification 
that we have performed is essential for one and two- 
dimensional Bose gases, it leads only to minor changes 
in the three-dimensional case. This can be seen by con- 
sidering the temperature dependence of the condensate 
density. At zero temperature, the fractional depletion 
that results from the Popov theory was first calculated 
by Lee and Yang |^ and equals {?> / "i) ^JnaFJ^i , where 
the T-matrix is taken to be r^^(— 2/i) = Airah^/m and 
a is the scattering length. The result that follows from 
Eqs. (|) and (|) is (32/3 - 2^/27r) ^/na^/n. The frac- 
tional depletion is thus approximately 2/3 of the Popov 
result. This is in fact the largest change in the conden- 
sate depletion, since the effects of the phase fluctuations 
decreases at larger temperatures. The critical tempera- 
ture Tbec is found by taking the limit no ^ in Eqs. (^) 
and (^. These expressions then reduce to the same ex- 
pressions for the density and chemical potential as in the 
Popov theory. This implies that our critical tempera- 
ture for Bose-Einstein condensation coincides with that 
obtained in the Popov theory and with that of an ideal 
Bose gas. 

Trapped Bose gases. — We now generalize our mean- 
field theory to a Bose gas in an external potential 
y™*(x). First, the right-hand side of Eq. is expressed 
in terms of the Bogoliubov coherence factors Uk and fk, 
which are then replaced by Uj(x) and Vjix). The latter 
satisfy the Bogoliubov-de Gennes equations P2,Pq] 



hujj- — V^- Ai(x)+ 2T2i^(-2^(x))n(x) 



2m 



Uj(x) 



+r2B(-2A*(x))no(x)w,(x) = , (10) 



huj, V^- ^(x) + 2T2J^(-2^(x))n(x) 

2m 



Wj(x) 



+r^"(-2^(x))no(x)w,(x) = , (11) 

where no(x) = |V'o(x)P and /i(x) — — V°^'^{tC}. Note 
that we have chosen Uj (x) and Vj (x) to be real. In some 
cases, e.g. when the macroscopic wavefunction "00 (x) 
contains a vortex, one cannot choose these amplitudes 
to be real and our equations must be generalized to in- 
corparate that fact. In terms of these particle and hole 
amplitudes, the expression for the total density in Eq. (Q) 
reads 

n(x) = no(x) + ^ S [wj(x) -l- Wj(x)]^ N{hujj) + Vj{x) 



r2B(_2^) = 



47rn^ 



1 



m \ii(2h /fima^ 



(9) 



x[uj(x) + Uj(x)] + 



r^^(-2M(x))no(x) 
2ej + 2/x(x) ' 



A,(x)]2 . (12) 
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Here, Cj are the eigenenergies and 0j(x) the eigenstates 
of the external trapping potential. In the large-j limit, 
we have that Uj{x) = 0j(x) and 



r^^(-2M(x))no(x) 
2e,- 



0,(x) . 



(13) 



As a result, the sum in the expression for the total density 
is ultraviolet finite since the second and third term cancel 
each other in the large-j limit. Furthermore, the inhomo- 
geneous generalization of Eq. (^) becomes the nonlinear 
Schrodinger equation 



Ai^o(x) 



2m 



y°'^'(x) 



r2B(_2A.(x))(2n'(x) + |^o(x)|2) 



V'o(x) . (14) 



Finally, the expression for the phase fluctuations in 
the case of a trapped Bose gas can be obtained from 
(x(x)x(x')), which is equal to the part of 



V"o(x)no(x') I L 



r2^(-2/.(xO)no(xO 
2ej + 2^(x') 



0,(x)<^,(x') , (15) 



which is symmetric under exchange of x and x'. Note 
that this correlation function is free of ultraviolet diver- 
gences in contrast to the long-wavelength result used in 
Refs. [011. 

Discussion. — We have proposed a new mean-field 
theory for dilute Bose gases in arbitrary dimensions, in 
which the phase fluctuations are treated exactly. We re- 
produce exact results in one dimension and the results in 
three dimensions are essentially the same as those pre- 
dicted by the Popov theory. The exact treatment of the 
phase fluctuations has solved the long-standing problem 
of infrared divergences in one and two dimensional Bose 
systems. This opens up the possibility to study in the 
same detail as in the three-dimensional case, the physics 
of low- dimensional atomic gases. As previously men- 
tioned, we can also incorporate the Kosterlitz-Thouless 
transition and for instance perform an ab initio calcu- 
lation of the critical temperature as a function of the 
scattering length and the density of the gas. 

The inhomogeneous generalization of the modified 
Popov theory describes in a selfconsistent manner the 
density profile and the full crossover from a condensate 
to quasicondensate in one and two-dimensional trapped 
Bose gases. In practice, however, we expect that a good 
approximation is obtained by simply calculating the 



densities no(x) and n'(x) in the Thomas-Fermi approxi- 
mation, i.e., by applying Eqs. (|^) and (||) locally at every 
point in space with a chemical potential equal to /i(x). 
Work in this direction is also in progress ||l9| . 
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